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Rasters Based on Hadamard Matrices for 
Oscillation Mode 

¥17F0163G Leningrad IZVESTIYA VYSSHIKH 
LCHEBNYKH ZAVEDENILY: PRIBOROSTROYENIY!I 
in Russian No 12, Dec 90 pp 74-77 


[Article by V.V. Afanasyev and V.B_ Shlishevskiy, 
Novosibirsk Institute of Geodetic, Aerial Photography, 
and Cartography Engineers] 

UDC 681,782.5 


[Abstract] Use of rasters based on Hadamard matrices 


for spectrometers operating with selective modulation of 


the luminous flux in the oscillation mode is considered, 
the essential requirement being that their autocorrela- 
tion function have no maxima besides those character- 
izing the dependence on the number of raster element. 
While such rasters constructed for selective modulation 
in the switching mode already have this property, their 
structure needs to be modified so as to make accurate 
quantitative measurements in the stable and reliable 
oscillation mode feasible. Such a modification involves 
transposition of rows and columns in the original matrix, 
which will not break its orthogonality. Sifting all possible 
transposition variants becomes simple when used 1s 
made of the relation H, = PH,Q, which represents the 
equivalence of two Hadamard matrices H, and H,. The 
monomial transposition matrices have each only one 
nonzero element in each row and in each column: a+! or 
a -l, representing respectively transparent and opaque 
elements of the raster. Figures 2; references 5. 


Measurement Data Processing Method and 
Equipment for Surface Shape Inspection 


917FO163F Leningrad IZVESTIYA VYSSHIKH 
UCHEBNYKH ZAVEDENILY 
PRIBOROSTROYENIYE No 12, Dec 90 pp 69-74 


{Article by M.A. Kropotkin (deceased), A.L. Melkonyan. 
and M.B. Stolbov, Leningrad Institute of Electrica! Engi- 
neering imen: V.I. Ulyanov (Lenin)] 


UDC 621.375 


[Abstract] One noteworthy application for optoelec- 
tronic systems is monitoring deformations of intricate 
surfaces in the product quality control process. Such a 
system consists essentially of an optoelectronic trans- 
ducer, a device for orienting the inspection object rela- 
tive to it, a scanner, a recording instrument, and an 
interface to a computer. Scanning and recording can be 
combined by using a photodiode matrix or a photoelec- 
tronic charge- coupled-device matrix. The operation of 
such a system is demonstrated on inspection of axisym- 
metric objects such as reflector antennas, paraboloids of 
revolution. The processing of measurement data 
depends on the description of the inspected surface in 
terms of approximating functions, Zernicke polynomials 
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having been selected for monitoring large deformations. 
The measurement procedure is accordingly determined 
by the parameters of the surface model which need to be 
estimated. The algorithm of calculating surface deforma- 
tions Pt points on the surtace takes into account the 
distance from scanner to inspected surtace. It includes 
correction of the object position relative to the trans- 
ducer, either for compensation of the positioning error 
or for minimization of the effect of surface deformations 
on the antenna pertormance characteristics. Such a cor- 
rection involves determination of th. approximating 
Surface, which is done Dy variation of the object position 
parameters B' = unordered set: 57,5x.dy dy*,dy” (linear 
displacements in x,y,and z directions, angular displace- 
ments by rotation about X and Y axes). The magnitudes 
of these variations are determined from the measure- 
ment data by the method of least squares and the effect 
of each of these displacements on the deviation of the 
object surface along the normal to it 1s described in terms 
of respective influencing functions. Compensation of 
homologous surtace deformations is considered, the nec- 
essary number of inspection points being established by 
expansion of the influencing functions into Zernicke 
polynomial series and compensation then being 
achieved by minimization of residual phase distortions 
in the antenna aperture. From the measurement data can 
also be determined the change in antenna performance 
parameters, gain and radiation pattern, owing to com- 
pensation of homologous surface distortions. Figures |; 
references 3 


Peculiarities of Reproduction Tomographic Images 
With Varifocal Mirror Objective 


91 7FOI63B Leningrad IZVESTIY.1 VYSSHIKH 
UCHEBNYKH ZAVEDENIY. PRIBOROSTROYENIYE 


in Russian No 12. Dec 90 pp §2-56 


[G.V. Mamchev. Novosibirsk Institute of Electrical 
Communications Engineering] 


UDC 621.397.13:681.3.04 


[Abstract] Reproduction three-dimensional multi- 
planar tomographic images in medical practice 1s con- 
sidered, use of a varifocal pancratic mirror objective 
having be 'n proposed for computer-aided tomography 
The process involves spatial superposition of two- 
dimensional images depicting each another parallel sec- 
tion of the object and to be seen one behind the other in 
Space in the following sequence: 1,2,...,n- 
I.nwn-l.n-1,....3,2,1, 2.....n.... (n- number of synthesized 
images of parallel sections) One possible version of a 
varifocal mirror objective 1s an aluminum-coated mylar 
film spread over a loudspeaker and splitting incident 
light by reflecting one part of it at a 90° angle onto the 
screen of the cathode-ray tube while transmitting the 
other part straight through. The design problem of 
determining that number n of images 1s tackled by taking 
into account the characteristics of binocular vision. The 
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depth of vision threshold needs to be determined accu- 
rately on a Statistical basis, a special experiment 
involving reception of televised images having for this 
purpose been performed in which the background lumi- 
nance was varied over the 2-40 cd/m? range and the 
signal-to-noise ratio was varied over the 6-40 dB range. 
The design procedure is demonstrated on such a mirror 
objective for a 61LK2B kinescope: aperture diameter of 
spherical mirror 2R = 610 mm equal to kinescope screen 
size (diagonal dimension), distance from kinescope to 
mirror in its central position a= 500 mm, inirror film 
vibrating at a frequency f, = 50 Hz with a maximum 
amplitude A,,, = 1-2 mm, n= 7-14, frequency of vertical 
scan fy = nf, = 350-700 Hz, number of resolution lines z 
= 625, frequency of line scan f, = 19-438 kHz. Figures 2; 
references 3. 


Processing Optical Measurements Data in 
Multiposition Small-Base Observation System 


917F0163D Leningrad IZVESTIYA VYSSHIKH 
UCHEBNYKH ZAVEDENIY: PRIBOROSTROYENIYE 
in Russian No 12, Dec 90 pp 60-65 


[Article by A.I. Daugavet and Ye.V. Postnikov, Lenin- 
grad Institute of Electrical Engineering imeni V.I. Uly- 
anov (Lenin)] 


UDC 621.391 


[Abstract] The performance of an optical multiposition 
observation system for determining the angular coordi- 
nates of remote objects is analyzed for accuracy of their 
estimates, that accuracy depending largely on the length 
of the ““measurement base’’. Systems with space diversity 
of at least two stations of the same type are considered 
for estimating the radius- vector of a point in space on 
the basis of the readings of its coordinates. The covariant 
estimation matrix is constructed ir such a system and 
the Rao-Kramer inequality is applied here, assuming a 
Mises distribution of the random quantity C, = (@,,@,;* 
on the [-1,1] segment (@,- readings of real direction 
angles w,", vectors w,- random points on a unit sphere) 
so that its probability density function is p(z) = Ae”. The 
results indicate that such a system must have a measure- 
ment base sufficiently long to make the angles between 
the directions from its stations to the target at least one 
order of magnitude larger than the angle measurement 
errors. References 3. 


Increasing Corrective Power of Spherochromatic 
Aberration Compensator and Its Use in 
High-Speed Photographic Cameras With 
Scanning by Mirror 

9/7F0163C Leningrad IZVESTIYA VYSSHIKH 
UCHEBNYKH ZAVEDENIY: PRIBOROSTROYENIYE 
in Russian No 12, Dec 90 pp 56-60 


[Article by A.V. Belinskiy, 1.V. Petrochenko, and A.V. 
Plokhov, Moscow Institute of Geodetic, Aerial Photog- 
raphy, and Cartography Engineers] 
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UDC 778,534.83 


[Abstract] Use of an afocal aberration compensator (L. 
Nefedov, OPTIKO- MEKHANICHESKAYA PROMY- 
SHLENNOST No 6, 1961) in high-speed photographic 
cameras in addition to the 30-200 lenses, depending on 
the number of frames to be recorded, is considered as a 
way to compensate the higher-order chromatic as well as 
spherical aberrations produced by switching of these 
lenses. Such a compensator must not only correct the 
angular errors of a lens but also suppress the secondary 
spectrum due to differences between the refractive 
indexes for different wavelengths. Following a perfor- 
mance analysis of a “thin’’ compensator for one lens, 
based on applicable formulas of geometrical optics, a 
design analysis of such a compensator yields dimen- 
sional and optical data necessary for selection of the 
compensator material. Such a compensator with a high 
corrective power in C, D, and F light has been designed 
with plain glass rather than difficult to machine KI- 
crystal and with a minimum of optical surfaces in the 
intermediate objective. Figures 2; tables 1; references 8. 


Mirror Compensator of Image Rotation 


917F0163E Leningrad IZVESTIYA VYSSHIKH 
UCHEBNYKH ZAVEDENIY: PRIBOROSTROYENIYE 
in Russian No 12, Dec 90 pp 65-68 


[Article by M.V. Dorofeyeva, State Institute of Applied 
Optics] 


UDC 535.317.2 


[Abstract] Considering that mirrors offer several advan- 
tages over prisms as compensators of image rotation in 
optical scanning systems, they are not aberrationless and 
less lossy but also easier to produce and at lower cost, 
their effectiveness when placed in the path of converging 
light rays between the objective and the plane of images 
is evaluated and compared with that of a prism compen- 
sator on the basis of their respective geometrical config- 
urations. Two factors are taken into account, that a 
prism cannot transmit a cone of converging rays with an 
angle wider than the critical angle @,., = tan (n/2k) (n- 
refractive index of prism material, k- ratio of the length 
of path of light rays in the prism to the diameter of a 
cylindrical light beam) and that the refractive index of a 
mirror compensator is n = |. A compensator consisting 
of three mirrors with a 120° angle between the first one 
and the third one has been designed for a cylindrical 
beam of light rays, with the same k = §.2 ratio as in an 
Abbe prism. A design and performance analysis of this 
compensator, its field of view depending on the aperture 
stop of the objective and the critical angle for a conical 
light beam being 15°, confirms its suitability and indi- 
cates indicates how it should be mounted. Figures 3; 
references 2. 
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Design of Mirror Objectives With Planoidal 
(Nearly Plane) Surfaces 


917F0163H Leningrad IZVESTIYA VYSSHIKH 
UCHEBNYKH ZAVEDENIY: PRIBOROSTROYENIYE 
in Russian No 12, Dec 90 pp 77-81 


[Article by G.I. Tsukanova, Leningrad Institute of Pre- 
cision Mechanics and Optics] 


UDC 535.833 


[Abstract] The design of three-mirror objectives with an 
almost aspherical and thus planoidal second or third 
mirror is analyzed, such a mirror having a large radius of 
curvature at the vertex and an infinitely large eccen- 
tricity squared. The analysis is based on the theory of 
third- order aberrations, these aberrations being usually 
compensated by surface deformation. Inasmuch as con- 
ventional three-mirror design methods are not appli- 
cable to planoidal mirrors, a new method is proposed 
which corrects not only spherical aberration but also 
coma and astigmatism. All three parallel mirrors, the 
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third mirror between the first and the second, are 
replaced with smooth surfaces and the following param- 
eters are stipulated: distances -d, from first mirror to 
second mirror and d, from third mirror to second 
mirror, heights h, > h, > h, of null ray on the principal 
planes of all three mirrors, distance -6 from from vertex 
of the second mirror to plane of the image, and angles a 
of null ray with the common optical axis. The null ray 
impinges on the first mirror parallel to the optical axis 
(a, = 0 and leaves the third mirror at angle a, = h,/f (f- 
second focal length), a, = -1 in objectives with no 
intermediate image or with two intermediate images and 
a, = | in objectives with one intermediate image after 
reflection by the first mirror or by two mirrors. The 
design procedure involves setting up and solving a 
system of equations which describes simultaneous full 
compensation of all three kinds of aberrations, in this 
case third-order aberrations of each kind. This system of 
equations yields expressions for the necessary deforma- 
tions of mirrors. Numerical calculations have been made 
for both meridional and sagittal sections in an objective 
with a 1:4 relative aperture and 3° field of view angle. 
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On Problem of Evaluating Wind Energy Potential 
of Armenia 


917FO183aYerevan IZVESTIYA AKADEMII NAUK 
ARMYANSKOY SSR: SERIYA TEKHNICHESKIKH 
NAUK in Russian No 4, 1990 pp 177-181 


[Article by A. A. Mardzhanyan, G. S. Petrosyan, and K. 
V. Khachatryan. First paragraph is authors’ abstract of 
article] 


UDC 621.311.24:621.548 


[Text] The specific energy characteristics of wind flow 
can be determined on the basis of Weibull distribution 
for regions promising from the viewpoint of developing 
wind energy. The parameters of the distribution can 
easily be calculated on the basis of statistical data on 
wind speeds using simple relations. A total of 1.5-2 
billion kW x hr of electric power per year can be 
generated if 1,500-2,000 modern wind generating plants 
are installed in promising regions of Armenia. 


Making the requirements more rigid on the ecological 
cleanliness of electric power generation and also the 
absence of local extractable types of fuel resources deter- 
mined the timeliness of developing alternative energy in 
Armenia in general and of wind energy in particular. 
Accurate analysis of the wind energy potential and 
development of a wind cadaster are a necessary base for 
development of wind energy in the republic. A prelimi- 
nary step of studies in this direction was developing a 
method of determining the energy characteristics of the 
wind flow at the wind vane level, developed at the 
Department of Electric Power Plants, Networks and 
Systems of YerPI [Yerevan Polytechnical Institute]. 


The basis of the method is approximation of the real 
probability density distribution of the wind speeds of a 
given locale using two-parameter Weibull distribution 





, k-! 2 
f(u) = = (+) exp -(=) |: (1) 
Cc c { c 


where u 1s wind speed at the height of the cadaster, m/s, 
k > Oandc > | are Weibull distribution parameters. 


4 detailed description of calculating the values of param- 
eters k and c on the basis of available statistical data on 
wind speeds for a given locale was presented in [1]. Let 
us dwell here on the problem of determining some 
specific energy characteristics of wind flow, which are of 
primary interest for solution of applied wind energy 
problems. 


One can show [2, 3] that the kinetic enjergy of the wind 
flow is proportional to the value u*, multiplied by the 
probability that wind with given velocity will occur 
Accordingly, to determine the value of wind speed with 
maximum energy, using expression (1), after simple 
transformations, we find 
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The value of u,,,. permits one to determine the maximum 
energy included per unit surface of wind flow in a given 
locale through the year [4] 


W,, = pu, f(u,,) 8760, (3) 


where f(u,,.) is the probability of wind velocity w.th 
maximum energy, determined according to (1), and p 1s 
the density of air. 


However, not all the energy determined by expression (3) 
can be used to generate electric power. In the ideal case, 
it is possible to use only 59.3 percent of the value of W,,,,. 
(the Betz limit [4]). Accordingly, the maximum electric 
power generated during the year with unit flow surface of 
an ideal wind turbine is determined according lo 


We =Com Wome Com = 0,593. (4) 


3 Pp 


The value of C,,, called the wind utilization factor, is less 
for real wind energy plants than the maximum value of 
C_.. and is equal to (0.3-0.47) as a function of the type 
of wind turbine [3, 4]. 
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Figure. Probability Density Distribution Curves of Wind 


Speed According to Weibull for Different Regions 


Key: |. Sisian Pass; 2. Pushkin Pass; 3. Yeratumber, 4 
Savan-ozernaya, 5. Aragats 


Based on available statistical information on the mean 
monthly values of wind speeds over a five-year period 
Weibuii distribution curves (figure) were plotted for five 
promising areas of the Armenian SSR, and the specific 
energy characteristics of wind flow were calculated 
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(table). As follows from the figure, the greatest shift of 
the maximum probabilities toward high wind speeds (> 
5 m/s) are observed for the Sisian Pass. The distribution 
curve for the Pushkin Pass is distinguished by a more 
sloping nature compared to the others and has the 


second shift of maximum probabilities toward the region 
of high velocities. The distribution curve for the Sevan- 
ozernaya region has a sharp peak for wind speeds of 4-5 
m/s. It follows from the table that the region of the Sisian 
Pass has the greatest wind energy potential. 





Table. Specific Energy Characteristics of Wind Flow 









































Observation point c, m/s Ume, m/s flume) Wom, kW x hr/m? Wam. kW x hr/m2. 
Sisian Pass 9.38/8.03 10.04/8.13 0.395/0.543 2,059/1,646 1 ,222/976 
Pushkin Pass 7.79/6.47 11.75/8.73 0.054/0.101 497/382 295/225 
Yeratumber 5.05/5.05 6.43/5.80 0.167/0.315 255/348 151/207 
Sevan-ozernaya 4.73/4.18 §.35/4.53 0.379/0.651 329/343 195/604 
Argats 4.11/4.03 5.95/5.46 0.120/0.160 143/148 85/88 














|. Calculations were made for normal atmospheric conditions. 

















2. Data for 1981 are given in the numerator, while mean data for 5 years are given in the denominator 





The reliability of calculations of the specific energy 
characteristics as a whole is determined by the accuracy 
of the initial meteorological information. The meteoro- 
logical stations supplying this information in the 
republic are located to assemble averaged indicators— 
they are frequently on the territory of populated points 
or near them, and are shaded by structures, local features 
of relief and so on. Moreover, the complex mountainous 
relief of Armenia is an extremely significant flow- 
forming factor, and a sharp deformation of the surface 
layer of wind occurs due to its influence. Air masses are 
held and deflected by mountain ridges, and the stream is 
constricted and as a result the flow is accelerated in the 
passes and troughs of the mdges. The wind speed can 
increase 2-2.5-fold at slight distance of several hundred 
meters, due to the characteristics of the relief. This local 
acceleration of the flow is very significant for generation 
of electric power, since the capacity of the flow increases 
in proportion to the cube of wind speed. Thus, regions 
with an abundant wind potential fall out of the visual 
field of the existing network of meteorological stations 
and a synchronous measuring system that permits one to 
determine the three-dimensional pattern of wind flow in 
the troughs of ridges and in passes must be developed for 
full evaluation of it. Moreover, the methods of taking 
into account micrometeorological features of the studied 
regions and of vanation of wind speed as a function of 
altitude must be worked out 


Having the specific energy characteristics of the wind 
flow of a given region, it 1s easy to calculate the annual 
generation of electric power, produced by using some of 
the existing industrial wind energy plants. Thus, let us 
consider a wind energy plant of the Nibe type (Denmark) 
with wind utilization factor of 0.42 and with air-swept 
surface of 1.257 m*. If a given wind energy plant is 
installed in the region of the Sisian Pass, the maximum 
annual generation of electric power will comprise | ,659 
MW x hr. Accordingly, if 600 of these wind energy plants 
are installed in this and other promising regions, annual 
generation of | billion kW x hr of electric power is 
possible. And if |,500-2,000 modern wind energy plants 





of medium capacities are located in promising regions of 
Armenia, | .5-2 billion kW x hr of electric power annually 
can be generated. This number of installed wind energy 
plants is not unrealistic if a state program is adopted for 
development of wind energy. Thus, 1,500 wind energy 
plants were installed in Denmark in 1986, and approxi- 
mately 200 new ones are turned over annually for 
operation, while more than 15,000 wind energy plants of 
medium capacities with total installed capacity of more 
than 600 MW were installed in the United States (Cali- 
fornia) for the same year. 
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The UkSSR Coal Industry in 1990 


917F0206A Kiev UGOL UKRAINY in Russian No 4 
1991 pp 55-62 


[Unattributed article under the “Chronicles” rubric 
“The UkSSR Coal Industry in 1990} 


[Text] In the past few years the republic's coal industry 
has experienced a reduction in its level of coal recovery 
and efficiency as a result of the aging of mine resources, 
an increase in the depth at which seams are developed. 
an increase in the complexity of geologic and mining 
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conditions, and unsatisfactory material and technical 
support of its mines. Various organizational flaws, a 
reduction in labor and industrial discipline, and a wors- 
ening of the social economic conditions of miners’ labor 
and daily living conditions have all had a negative effect 
on work in the sector. The system of measures stipulated 
hy government decree No. 608 has not been fully imple- 
mented. This has created tension in the republic’s coal 


JPRS-UEQ-91-008 
27 June 1991 


basins. All of these shortcomings reduced the sector's 
operating efficiency in 1990. 


Compared with 1989, the recovery of all coals in 1990 has 
been reduced by 4.47 million tons, amounting to 175.7 
million tons. This in turn is only 93.8% of the planned figure 
(Table 1). Except for the Aleksandriyaugol, no association 
has fulfilled its coal recovery plan. Mos? are noticeably 
below their yearly coal recovery volume. 


































































































7 Table 1. 
Production Recovery of Ati Coals | Recovery of Coals for Coking 
Association | | a 
Pianned, | Actual, Percentage | Percentage | Planned, | Actual, Percentage Percentage 
thousands of | thousands of _ of Plan | of 1989 | thousands of _ thousands of of Plan of 1989 
7 ; | tons | tons | it Amount | tons | tons | Amount _ 
Donetsk Oblast | 85,148 | = 80,451 99.5 | 91.7 | 49.899 | 47,754 | 95.7 91 
Donetskugol | 19,160 | 18.728 | 97.7 «| 927 14,572 | 14301 98 1 90.7 
Makeyevugol = | 11,930 | 11.207, «|| 93.9 | 8 11,365 | 10.733 | = 944 91.7 
Krasnoarmeyskugol | 8.690 | 8.462 | 97.4 | 89.8 | 8,300 | 8.079 if 97.3 | 90.3 7 
Selidovugol | 315 | 5.184 | 97.5 | 93.2 f | 
—_——_ + + + + + 7 + ————_—_—_____+ ——___—__— 
Dobropolyeugol 6,955 6,817 | 9.0 | 97.1 | 5,440 5,328 97.9 99 
Asiomvdie | $963 | 5390 | 04 | 852. =| $412 | 4893 | 1904 | 47 
—— _ + . + + + — 4+—_______ —_—+___--___ —____—_ 
Dzerzhinskugol | 2,840 | 2,591 | 91.2 | 88.2 | 2,84( | 2,591 9] 2 it 88 2 
Ordzhonikidzeugol | 4555 | 4.086 | 89.7 | 85.) | 1.570 | ] 433 | 9! 2 _| 91.7 7 
Shakhterskugol | 4870 | 4560 | 936 | 971 | | 
Oktyabrugo! | 6,630 | 5.987 | 90.3 | 93.2 | 40K | 3964 | 9 | | 9? 7 
Torezantratsit 8,240 If 7.439 | 90.3 | 91.7 | | | 
Lugansk Oblat | 55.975 | 51619 | 922 | #«2©931 +«+| 12.4890 | 11.813 949 99 7 
Luganskugol | 8.400 8 004 il 9§ 3 ia 93.7 | 400 | 49 ~ 9 T 4) 
Stakhanovugol 6305 | 6017 | 954 | 987 | 3,005 | 287 T0955 | 98s 
Pervomayskugol | 3.850 | 3,242 | 842 | 924 «| ~~ 1,595 1307 |) 3819)ts« 
Lisichanskugol 3280 | 3162 | 4964 +| «+1029 | if Ty | 
Krasnodorugol | 7,600. «| ~——~*7,362 %9 | 1002 | 7.45 1267 | #975 «| «10 
D nbessantreteit | 7,000 iz 6,080 | 86.9 | 85.4 | | 
Antratsit. = 4140 | 3,432 29 | Bi | | — 
Rovenkiantratsit. | 7,200 | ~=6.773. | sO 18 | | - | 
Sverdlovantratsit | 8,200 7.54) T 92 | 92.2 ~T | 
Pav logradugol i3sis | 1298 | 961 | #42»908 | +3003 | 3102 | 03 
Ukrzapadugol | —s11,830 «| ~=SO«10.824 | 2Sti90SC*«SYté‘<~CSD*SS*~‘St‘<i«‘“SSCIN~C*dSS*CICSt*~<‘ia SSSCS*«é SS 
\ieksandnyaugol | 9.230 «|~S=«OSDS].S«)SSC«*O a6. TT. | 
Thr ughout | 175.698 164.810 T 938 91.5 | 71,162 | 67,787 | 9S | 
UkSSR | | | iL | | 





T 


ihe underground method was used to recover 158.3 
million tons of coal (93.3% of the plan and 91.3% of the 
|989 level). and the open-pit method was used to recover 
“5 millon tons of coal (106.5 and 97.1%, respectively) 
The hydraulic method was used to recover 2.8 million 
tons (105 and 98%). The average daily recovery of coal in 
the republic equaled 480,000 tons as opposed to the 
figure of $18,000 tons stipulated in the plan. Recovery 
n | square meter of the web area of a seam amounted 


to 1.78 tons, which is 7.5% above the indicator stipulated 
in the plan and 0.6% of the 1981 indicator 


Coals for coking were recovered in the amount of 67.8 
million tons or 95.3% of the amount stipulated in the 
plan and 93.3% of the 1989 level. The only association to 
meet the plan for recovery of coal for coking was the 


Paviogradugol Association, which 1s located in an area 
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where recovery of such coals has increased in recent 
years. It has, however, dropped noticeably in the central 
region of the Donbass at the Pervomayskugol and Uk:za- 
padugol associations 


In 1990, 206.5 million tons of coal was delivered to con- 
sumers. This was 2.2 million tons short of the plan. Man\ 
mines and associations did not meet their contractual 
obligations. The ash content of the coal recovered amounted 
to 29.6% with a standard of 29.6%. The ash content of coal 
delivered was 17.8% versus a norm of 18.4% and was 0.6% 
lower than in 1989. Balances at warehouses holding all coals 


equaled 4,720,000 tons as opposed to the standard-specified 
amount of 4,240,000 tons. The balances of coals for coking 
amounted to 846.000 tons versus the 1,090,000 tons called 


for in the standard 
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Prepared coal reserves in sections amounted to 
1,947,000 tons, whereas the plan called for 2,115,000 
tons; 785,000 tons was prepared for removal, whereas 
the plan called for preparing 900,000 tons for removal 
The amount of stripping operations was lower than in 
1989. 


Table 2 gives performance indicators both for all devel- 
opment openings (by the fiscal and contractual methods) 
and for stripping and development operations. The 
figure for total developments fell 3% short of the plan 
and 3.7% short of the 1989 level. The corresponding 
figures for stripping and development operations were 
7.1 and 9.2%, respectively. 





Table 2. 





Production 


Development Openings (by the fiscal and contractual methods) 





Association 


Including Stripping and Development 
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Ordzhonikidzeugol, Shakhterskugol, and Aleksandriyaugol) 
met their plan indicators for both development openings 
and stripping and development workings 


plan for opening developments met at the Selidovugol, 
Dobropolyeugol, Oktyabrugol, and Lisichanskugol associa- 
tions were especially low. Only three associations (the 































































































Table 3. 
Production Association | Average Active Number of Average Monthly of Face Average Active Line of Average Daily Load on an 
Cleared Faces Advance, m | Cleared Face. km Active Cleared Face, t 
7 T am a 
Donetsk Oblast | 894.4 28.6 130.4 263 
Donetskugol 191.5 26.4 32.2 295 
T - T T — 
Makeyevugol 121.9 | 25.9 19.5 265 
Krasnoarmeyskugol 33.7 $0.2 62.6 766 
Selidovugol 33.0 46.1 49 464 
Dobropolyeugol | 33.2 66.8 46.1 676 
Artemugol 132.9 23.1 14.3 2 
Dzerzhinskugol | 68.5 19.5 79 7 
Ordzhonikidzeugol 105.6 22.8 11.3 7 
Shakhterskugo! 41.8 25.3 8.2 285 
Okty abrugol §4.0 28.7 10.1 33¢ 
Torezantratsit 78.3 25.8 120 27 
Lugansk Oblast 459 28.7 77.6 338 
——— He —- 7 + —— — 
| uganskugol 66.4 33.6 10.9 3 $¢ 
Stakhanovugo! 95.7 20.3 15.2 18 
— 7 a + + + + + —--— 
Pervomayskug 37.4 21.1 68 272 
~ - - +— — +> — — —_——_—— — — 
l hanskug 22.4 33.5 39 473 
Krasnodonug $6.7 | 33.4 9.0 386 
- a + + —+— 7s — = 
Donbassantratsit 67.6 240 11.6 295 
—— —_+- + +— + —_—— 
Antrats 31.7 21.3 §.7 342 
Rovenkiantratsit 40.4 35.9 7.6 495 
sonhtibca 1 in 4 sl pieeeasiepotennioneae 
Sverdiovantratsit 40.8 48 3 6.7 $72 
Pav logradugol 69.5 62.9 11.8 S68 
———__—__—- + + +> > - - 
L krzapadugol 749 417 111 4¢ 
- — +> 4 —+ - ——+ --- - . 
Aleksandnyaugol 16.5 417 14 SRI 
- - + — + —— 
Throughout UkSSR 1,409.8 | 30.9 226.3 31) 





Table 3 presents technical mining indicators with respect 
to all cleared faces and to working faces equipped with 
mechanized systems (abbreviated KMZ). Compared 
with 1989, the average active number of cleared faces in 
1990 decreased by 53. Their average active line was 
reduced by 5.3 km, and the average monthly face 
advance was reduced by 2.4 m. The average daily 
recovery from an active cleared face was less than 8 tons. 
The average active number of working faces equipped 
with mechanized systems remained practically at the 
1989 level, their average active line increased by 1.9 km. 
The average daily load was 1.6% above the amount 
specified in the plan but 3.4% less than the 1989 level 


In 1990 the relative share of workings made with mech- 


anized loading of coal and ore amounted to the following 
(the data presented in parentheses include cutter- 
loaders): throughout the UkSSR, 81.3% (33%); through 
the Donetsk Oblast, 83.9% (32.1%); for the Donetskugo!, 


85.8% (37.4%): for the Makeyevugol, 76.8% (30.2%); for 
the Krasnoarmeyskugol, 98.5% (79%): for the Sell- 
dovugol, 80.6% (43.6%); for the Dobropolyeugol, 96.7% 
(74%): for the Artemugol, 100% (0.9%): for the Dzerzhin- 
skugol, 99.6% (0.6%); for the Ordzhonikidzeugol, 92.1% 
(0%); for the Shakhterskugol, 59.0% (23.4%): for the 
Oktyabrugol, 64.1% (16.6%); for the Torezantratsit 
67.4% (6.8%); throughout the Lugansk Oblast, 72.2% 
(12%); for the Luganskugol, 80.4% (20.1%): for the 
Stakhanovugol, 76.5% (3.4%); for the Pervomayskugol. 
66.0% (10.3%); for the Lichanskugol, 86.1% (37.3%): for 
the Krasnodonugol, 83.4% (29.2%): for the Donbassan- 
tratsit, 81.8% (11.6%): for the Antratsit, 59.0% (1.9%): 
for the Rovenkiantratsit 48.1% (6.3%); for the Sverdlo- 
vantratsit, 65.0% (1.2%). for the Pavlogradugol, 98.0% 
(96.8%); for the Ukrzapadugol, 86.5% (68.9%): and for 
the Aleksandriyaugol, 74.3% (74.3%) 


Table 4 presents key operating indicators for the year for the 


republic's coal enrichment facilities. A total of 132.02 
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tf coal was processed, and 81.65 million tons of 
concentrate was obtained. The yield of coals belonging to 
the coarse and medium classes amounted to 19.66 million 
tons, including 12.81 tons of anthracites. About 4 million 


tons cf coal briquettes were produced. As is evident from 


million tons ¢ 
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Table 4, however, the figures stipulated in the plan were not 
reached. Furthermore, a decrease in all indicators when 
compared with the data for 1989 was noted. This was largely 
due io a decrease in the volumes of coal that was recovered 
and sent for enrichment. 















































Table 4. 

Indicator Enrichment of coals OO 

Planned, thousands of Actual, thousands of | Percentage of Plan | Percentage of 1989 
ee CS tons tons | ___ Amount | 
Coal processing hment plants 42,076 132,023 H 92.9 mi 94.1 
Including for coking 66.772 $4,340 | 95.7 T 740 
Dutput of rn nr ery 84.646 | 94.7 | 94.6 
Including for coking 3$,§7 34,544 | 97. 97.7. 
O 7 7 ore anc Medium-Class ¢ als —_ 21,162 19.661 | 92.9 O11 
including anthracites 14,001 12.81 | 91.5 91.¢ 
Coal processing at mechanized spoil-picking | 11,306 11,061 | 97.8 | 944 
nstaiiat S 
Production of coal bnquenes~S~S*«dYSt*«~<~«SD Ts a 





Enterprises’ production of consumer products at optimum 
prices was 4.3% below the amount stipulated in the plan. It 
was 4.6% lower in calculated prices. The production cost of 
| ton coal in January-November 1990 amounted to 34.87 
rubles versus the 32.47 rubles stipulated in the plan, and it 
was 5.94 rubles higher than the amount for the same period 
yf 1989. The average monthly number of industnal and 
production personnel in the republic’s coal industry was 
15,600 less than the number specified in the plan and 5,500 
less than the 1989 figure. The output of consumer product 
per worker was as follows: throughout the entire UkSSR. 
8.551 rubles versus the the 8,758 rubles stipulated in the 
plan; throughout the Donetsk Oblast, 7,445 rubles (7,600 
rubles in the plan); and throughout the Lugansk Oblast 
9.219 rubles (9,506 rubles in the plan) 


In 1990 production facilities to recover coal in the amount 
of 4.870.000 tons per year were introduced into the coal 


industry of the UkSSR. The first phase of the Krasnoarm- 
eyskaya-Zapadnaya mine of the Krasnoarmeyskugo! Pro- 
duction Association, which has a capacity of 1,500,000 vy 
was put into operation. Also put into operation were the 
second phases of the following mines: the Shakhterskaya- 
Glubokaya (Shakhterskugol), which has a capacity of 
1,200,000 t/y, the Komsomolets Donbassa (Oktyabrugol) 
which has a capacity of 550,000 t/y; and the Luganskaya No 
| (Luganskugol), which has a capacity of 500,000 vy. As a 
result of reconstruction, the capacity of the mine Komso- 
molskaya (Antratsit) increased by 400,000 tons. and that of 
the mine Krasnyy partizan (Sverdlovantratsit) increased by 
100,000 tons. The Konstantinovskiy section of the associa- 
tion Aleksandriyaugol (capacity, 500,000 tons of coal per 
year) was also put into operation 


COPYRIGHT: Ukrainskoye respublikanskoye uprov- 
leniye VNTO Gornoye, 1991 
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Dynamics of Asymmetric Interaction Involving 
Body of Revolution and Rigid Barrier 
9/7F0196A Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 31-35 


[Article by I1.Ye. Khorev and N.T. Yugov, Scientific 
Research Institute of Applied Mathematics and 
Mechanics, Tomsk] 


UDC 539.3 


[Abstract] Asymmetric interaction of an undeformable 
barrier and steel cylinders striking it obliquely is consid- 
ered, of concern being the dynamics of this interaction 
and specifically the energy transfer during impact. The 
interaction is analyzed on the basis of numerical calcu- 
lations made by the method of finite elements in accor 
dance with the elastoplastic model of continuous media. 
Such calculations were made for three cylinders with the 
same initial diameter (D = 0.0125 m) but different 
length-to-diameter ratios (L/D = |, 2, 3). Each cylinder 
was made to strike the barrier obliquely at various angles 
(15°, 30°, 45°, 60°, 75°) and with various velocities (V = 
100 m/s, 300 m/s, 500 m/s) at each angle. The analysis 
covers a transient period o: 200 us. It reveals how the 
deformation of a cylinder as well as its pressure force and 
loss of kinetic energy as functions of time during that 
period depend not only on its initial relative length and 
its velocity but also on the angle of impact. The results of 
this analysis give a better insight into the mechanism of 
ricochet. Figures 4; tables 2; references 10. 


Experimental Study of Dynamic Deformation of 
Hollow Elastic Cylindrical Shells Immersed in 
Water and Vacuumized 


917F0196B Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 49-55 


[Article by I.I]. Anikyev, M.I. Mikhaylova, and Ye.A. 
Sushchenko, Institute of Mechanics, UkSSR Academy of 
Sciences, Kiev] 


UDC 539.3 


[Abstract] An experimental study of thin hollow circular 
elastic cylindrical shells under static pressure produced 
by immersion in water and vacuumization was made, of 
interest being their dynamic deformation by a weak and 
short subcritical pressure pulse under these conditions. 
The experiment was performed with 330 mm long and 
120 mm in diameter, 0.4 mm thick shells of DI6T 
duralumin turned on a lathe and 0.8 mm thick shells of 
glass-plastic (seven layers of E-3-10 glass cloth tape with 
EF 32-301-16 binder: ED-16 epoxy resin + bakelite 
varnish + acetone). Their walls were sealed by coating of 
shellac on the outside surface and their both ends were 
tightly closed by plugs, a pipe connected to a suction 
pump passing through one of them. After having been 
tested in air for residual strains, each sheil was immersed 
to certain depth within the 1|.1-1.5 m range in water 
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filling a 3 m high cylindrical tank 3.2 m in diameter and 
air was from its cavity. A pressure wave in water was 
then generated by explosion of a 15 mm long copper wire 
0.1 mm in diameter, through which a 48 pF - 2.5 kV 
battery of IKM-25-12-U4 noninductive capacitors had 
been discharged. The amplitude of the acoustic wave was 
about 0.3 MPa, varying from test to test but only within 
+/-3.5 %, and its steep front at a distance of 0.5 mm from 
the wire was followed by an exponential pressure fall 
with a time constant of about 50 ys. For measurements 
were used KFSP1-1-200A-12 strain gages with a | mm 
base or KF5P4-3-100V-12 strain gages with a 3 mm base 
in a TSTM-5 digital tensometer bridge. The error of 
pressure settings did not exceed 5 % and the error of 
strain readings was accurate within about 10 %. Pressure 
and strains were were also recorded on an S8-1 3 oscillo- 
graph over a transient period of 1.2 ms. All shells were 
found to behave as imperfect ones, hydrostatic pressure 
thus produced in them flexural as well as membrane 
strains. No strain of any shell exceeded 0.003 within the 
elastic range. The results indicate that a weak and short 
pressure pulse causes the dynamic membrane strains of 
such shells under a static load to alternate, their ampli- 
tude being larger and their duration being longer in the 
first compression half-cycle than in the following ones, 
and the magnitude of dynamic flexural strains during the 
tension-compression cycle depends largely on the total 
Static pressure load. Figures 4; references 10. 


Nonlinear Deformation of Cylindrical Shell by 
Local Thermal Shock 

9/7F0196C Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 55-60 


[Article by N.I. Obodan and N.B. Makarenko, Dnepro- 
petrovsk University] 


UDC 539.3 


[Abstract] Nonaxisymmetric deformation of thin-walled 
cylindrical shells and panels by fast local heating under 
uniform external pressure is analyzed on the basis a 
geometrically and physically nonlinear model which 
combines the theory of flat shells with the incremental 
deformation theory of plasticity. The total strain is 
accordingly expressed as the algebraic sum of its elastic, 
thermal, and plastic components. The two equations, of 
motion and of compatibility of strains, are solved for 
homogeneous initial conditions and boundary condi- 
tions of circumferential periodicity (closed shells) or 
boundary conditions at the edges (cylindrical par |s). 
The transient temperature field in a shell is found by 
solving Fourier’s parabolic second-order partial differen- 
tial equation of heat conduction for zero initial condi- 
tions and boundary conditions of convective heat 
transfer to the ambient medium, from a bare lateral 
surface and insulated bases. Integration of this equation 
yields two components of the temperature field: temper- 
ature of the median surface and temperature difference 
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across the wall thickness. A system of ordinary differen- 
tial equations is then constructed for deflections and 
both components of its temperature field. This system of 
equations being can be solved numerically by the Runge- 
Kutta method, after reduction to the Cauchy problem 
and determination of the initial displacement vector for 
each time step by Newton's method. The problem was 
solved for cylindrical steel shells 0.1 mm thick and 2.0 
cm long with a 4.3 cm radius first under zero static 
pressure, then successively under static pressures equal 
to 60 % and 85 % critical, while being heated by a locally 
incident thermal flux at characteristic rates of 5x10" °C/s 
(n = 0,2,4,6,8). The results indicate that plastic deforma- 
tion begins sooner as the heating rate is increased, the 
region of plastic deformation expanding slowly when the 
heating rate is low (quasi-static deformation) and 
spreading instantaneously over the entire thickness of 
the shell or panel wall when the heating rate is high 
(dynamic deformation). In the latter case the size of the 
heat spot and thus the degree of thermal flux localization 
is not an influencing factor. Intermediate heating rates 
thus cover the range associated with strong interplay of 
comparable static and dynamic deformation effects. 
Figures 5; references 12 


Transient Processes in Electrically Excited 
Multilayer Piezoceramic Stack 


917F0196D Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 61-66 


[Article by A.E. Babayev and V.G. Savin, Institute of 
Mechanics, UkSSR Academy of Sciences, Kiev] 


UDC §39.6.013.42 


[Abstract] Transient processes in a piezoelectric stack of 
N plane layers following its excitation by an electric 
pulse signal are analyzed on the basis of the linear theory 
of electroelasticity. All layers have the same thickness h 
and identic:! mechanical properties. The lower surface 
of the bottom layer adheres rigidly to a base (z = 0) and 
the upper surface of the top layer (z = Nh) ts load-free, 
but the electric pulse impinges on that surface of the top 
layer (k = N) and then continues passing through succes- 
sive interlayer boundaries (layers k = N-1,...,2,1). The 
corresponding three systems of k equations (k = 

,2,...,.N), two equations involving normal displacements 
and electric potentials and one equation involving 
normal stresses and electric potentials for each layer, are 
solved for boundary conditions of unbroken contact 
between layers, zero normal displacement of lower sur- 
face of bottom layer (k = |), compatible normal displace- 
ments of all interlayer boundaries and compatible 
normal stresses at all interlayer boundaries (k = 1,2,...,N- 
1), zero normal stress at the upper surface of top layer (k 
= N), electric potentials ¥"’\. v,.., = (-1)*'Q(t)H(t) 
and “{"2 = & = (~1)* Q(t)H(t) (k = 1,2,....N). Here function 
Q(t) characterizes the waveform of the incident electric 
signal and H(t) is Heaviside’s unit function. Following 
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elimination of the electric potential from the two equa- 
tions where it appears with displacements, the remaining 
two systems of equations are solved by the Laplace 
transformation method. The results of a numerical anal- 
ysis for a piezoelectric transducer made of TsTBS-3 (B, 
Sr zirconate-titanate) indicate how the variation of its 
stress and strains depends on the number of layers and 
the signal waveform. Figures 5, references 5 


Nonaxisymmetric Stressed and Strained State of 
Anisotropic Multilayer Shells of Revolution 
917FO196E Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 11, Nov 90 pp 66-70 


[Article by G.M. Kulikov, Tambov Institute of Chemical 
Apparatus Design] 


UDC 539.3 


[Abstract] A stress and strain analysis is performed for 
nonorthotropic shells of revolution with uniform thick- 
ness consisting of N anisotropic layers, each layer having 
a different uniform thickness. The reference surface is 
placed in a system of curvilinear orthogonal coordinates 
a, and a,, its curvature and the Lame parameter A, then 
being functions of only one variable a, . The resolvent 
system of ten partial differential equations is supple- 
mented with four differentia! and four algebraic rela- 
tions for the relevant shell constants and variables, also 
with a fourth-order matrix equation for forces and 
moments. The system is then solved for zero boundary 
conditions identical to the classical boundary conditions 
in the Kirchhoff-Love theory of shells. The solution 1s 
expressed in the form of two vectors X(a, and a,), both 
as well as the load q assumed to be expandable into 
complete Fourier series on the [-x,n] interval. That 
system of partial differential equations thus reduces to a 
one of ordinary differential ones which can be integrated 
by the method of discrete orthogonalization, those four 
differential relations being reduced to algebraic ones 
The algorithm of calculating the deflection in this way 
has been programmed in the PL/1(0) version of the 
high-level PL/1 language in the form of the main ANST3 
procedure, for execution on a Standard System 1060 
computer. The procedure was tested first on a homoge- 
neous transversely isotropic cylindrical shell with rigid 
constraints at both ends while under a normal load q = 
q,cosna,, then on a cylindrical shell consisting of two 
orthogonally wound boron-plastic layers with similar 
constraints at both ends while under a similar load 
Figures 2; tables |; references 7 


Dispersion Equations for Electroelastic Shear 
Waves in Multilayer Media With Periodic 
Structure 

917FO196F Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 84-93 


[Article by L.P. Zinchuk, A.N. Podlipenets, and N.A 
Shulga, Institute of Mechanics, UkSSR Academy of 
Sciences, Kiev] 
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UDC 539,3:534. 1 


[Abstract] Dispersion equat ons for surface acoustoelec- 
tric shear waves in multilayer media with periodic struc- 
tures are derived, media considered here being those in 
the crystallographic 6mm class such as a wurtzite crystal 
or polarized piezoceramics. Propagation of such a wave 
through such a medium in an xy plane 1s described by a 
system of six partial differential equations for the stress 
tensor components o,, and a,,, the displacement w, the 
electric-displacement vector components D, and D,, 
and the electrostatic potential ». The solution to it 1s 
sought in the form unordered set: D, (x,y,t), 6,,(x,y,t), 
p(x,y,t), w(x,y,t) = unordered set: q, (y), Gly), p,(y), 
p(y) e'**'@" (t- time, y- radian frequency of wave, k- 
wave number). For this the system of equations 1s 
transformed into one of ordinary differential ones for 
vectors q = col(q,, q, and p=col(p,, p,. The monodromy 
matrix U(h) of this h-periodic Hamiltonian system has a 
characteristic equation det{U(h)- KE] = O which, 
according to the Lyapunov-Poincare recurrence the- 
orem, is K* + a,K’ + a,K° + a,k + | = 0. The dispersion 
equations for surface shear waves in such a medium are 
constructed on this basis, for a medium occupying the 
upper half-space y > 0 with a,,,,,.,, and @(x,0,t) = 0 
conditions at its boundary y = 0 with vacuum below. For 
a dispersion analysis of waves in media with specific 
structures, one can evaluate both the monodromy matrix 
U(h) and its matritzant U(y) in that system of ordinary 
differential equations or solve the Cauchy problem for 
that system numerically by the Runge-Kutta method. 
Such an analysis 1s demonstrated on three examples: |) 
surface waves in a structure formed by alternating layers 
of TsTC-19 (Sr zirconate-titanate) piezoceramic and 
epoxy resin, 2) surface waves in a structure formed by 
alternating layers of crystalline CdS and 7nO, 3) normal 
waves in a structure formed by layers of crystalline CdS 
and ZnO. Figures 6; tables |; references 23. 


Use of Asymptotic Method for Analysis of 
Finite-Amplitude Vibrations of Flat Shells 
917F0196G Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 93-99 


[Article by N.I. Zhinzher and V.Ye. Khromatov, 
Moscow Institute of Energetics] 


UDC 5$39.3:534.1 


[Abstract] Nonlinear free vibrations of thin elastic and 
isotropic shells, rectangular in the plan view with sides 
a,,a, and thickness h, are described by Andrianov’s 
simplified integrodifferential equations (1. V. Andrianov, 
PRIKLADNAYA MATEMATIKA I MEKHANIKA 
Vol 50 No |, 1986) analogous to Berger's equations for 
plates but with the median surface referred to system of 
curvilinear coordinates x,,x, coinciding with the prin- 
cipal curvature lines. This system of equations 1s solved 
by Bolotin’s asymptotic method (V.V. Bolotin, 
“Random Vibrations of Elastic Systems”, Izd. “Nauka” 
1979: IZV. AKADEMII NAUK SSSR: MEKHANIKA I 
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MASHINOSTROYENIYE No 4, 1962) for small vibra- 
tions far from the edges of the shell and thus non- 
momental ones. Insertion of the generating solution for 
the deflection function within this region into that 
system system of equations yields a resolvent equation 
which relates the frequency of nonlinear vibrations io 
their amplitudes and wave number through elliptic func- 
tions. Numerical calculations by this method reveal how 
the distribution of natural frequencies depends on the 
amplitude of vibrations, as demonstrated on a plate. 
Figures 5; references 8. 


Forced oy Oscillations of Long Pipe 
Containing Liquid 

917F0196H Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 11, Nov 90 pp 114-118 


[Article by V.S. Tikhonov. V.I. Komissarenko, and B.V. 
Yemets, Central Scientific Research Institute of Geolog- 
ical Exploration, Moscow] 


UDC 534.11:622.648 


[Abstract] Equations of motion are derived to describe 
longitudinal oscillations of a long pipe driven into the 
ocean bed for exploration, such a pipe containing a 
column of liquid above the mud which plugs it at the 
bottom. The equations are based on a linear relation 
between the change of pressure and the bulk train of the 
liquid column, assuming that the liquid is a linearly 
compressible one, and an analogous relation for the 
linearly elastic pipe in accordance with Hooke’s law. To 
the two ordinary differential equations of statics for 
infinitesimally small pipe and liquid column segments 
are added two partial differential equations of dynamics 
for the longitudinal forces acting on them in accordance 
with Newtons’s second law, including not only the pres- 
sure force and gravity but also viscous friction against 
surrounding sea water. Transformation of these equa- 
tions in four real variables (x- displacement of pipe, X- 
displacement of liquid column, AF- increment of pres- 
sure force on liquid, AT- increment of tension force on 
pipe) as functions of time into equations in corre- 
sponding four complex variables as functions of fre- 
quency reduces them to a linear boundary-value problem 
for a system of four ordinary differential equations 
which can be solved numerically by the ranging method 
This has been done for a 6000 m long steel drill column 
with a 0.146 m inside diameter and a 0.168 m outside 
diameter, weighing 41.5 kg/m and containing water, 
assuming a w(W) = const amplitude-frequency character- 
istic and a w(w) = O phase-frequency characteristic of 
platform heaving. Figures |; references 4. 


Propagation of Shear Waves Through Nonlinearly 
Elastic Body 

917F016401 Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 1, Jan 91 pp 127-129 


[Article by V.I. Yerofeyev and I.G. Raskin, Institute of 
Machine Science, USSR Academy of Sciences, Nizhne- 
gorodskiy branch] 
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UDC 539.3 


[Abstract] Propagation oi plane shear waves through an 
isotropic nonlinearly elastic body is described by a wave 
equation with a coefficient a in the a/p(Su/5x)?8*u/6x" 
term (u- transverse displacement) which additively com- 
bines the effects of quadratic, cubic, and quartic elas- 
ticity nonlinearities. The solution to this equation is, in 
the linear approximation, represented as a series of 
harmonics whose frequencies @, and wave numbers k, 
are related through the dispersion equation @,- c, k, (c, = 
(u/p)? denoting the velocity of shear waves, p- density of 
material, p- Lame constant of second order). In the case 
of a strong cubic nonlinearity the effect of self-action 
dominates over generation of higher-order harmonics so 
that the latter may be ignored. The solution can in this 
case be sought in the form of a single harmonic with a 
complex amplitude which slowly varies in space and in 
time. It is found by averaging over “fast” variables. A 
random process in a body with cubic nonlinearity is 
considered, a Gaussian noise for specificity, and its time 
correlation function B(t,x) is found in a form which 
facilitates solution of the following inverse problem: to 
calculate the coefficient of cubic nonlinearity from the 
change in spectral width of a noisy shear wave. Numer- 
ical calculations were made for a noisy shear wave with 
a frequency of about 5 MHz propagating through a 
specimen of cubically nonlinear D16T aluminum alloy 
with a strain 6u/5x =~ 0.0003 and a Aw,/Aw, = 1.1 or 10 
% change in the width of its spectrum, the nonlinearity 
coefficient being a = 10’ N/mm?. The authors thank A.1. 
Vesnitskiy for discussing the results. References 5. 


Longitudinal Damped Vibrations of Plates 


917F0164H Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 1, Jan 91 pp 124-127 


[Article by A.S. Kosmodamianskiy, Donetsk University, 
and M.N. Gofman, Mariupol Institute of Metallurgy] 


UDC $39.3 


[Abstract] Longitudinal vibrations of a plate with arbi- 
trary smooth outer and inner contours moving under a 
periodically fluctuating pressure while overcoming a 
resistance force proportional to its velocity are analyzed 
by Kosmodamianskiy’s approximate method for steady- 
state vibrations of multiply-connected bodies (A.S. Kos- 
modamianskiy, VESTNIK UkSSR AKADEMII NAUK 
No 4, 1988). Inasmuch as damped vibrations are of 
concern here, the two components of the displacement 
vector U(x,y,t), V(x,y,t) in the system of two partial 
differential equations of motion relating gradients of 
Stresses GO, ,, t,,) to components of vibration velocity 
adU/d5t, adV/5t and components of vibration accelera- 
tion pd,U/t*, p&°V/8t? are expressed in the form 
A(x.y)e"" ~ 2" (@- frequency of pressure fluctuations, a- 
damping coefficient). After insertion of Hooke’s law and 
introduction of the complex variable U )longitudinal 
vibrations, this system of equations is solved approxi- 
mately by that method for boundary conditions 
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depending on the problem. Both the general solution 
U°,, and the particular solution U*, (2,2) are obtained 
the latter in the form of a double integral. The damping 
of vibrations, in terms of displacements as functions of 
time, has been analyzed numerically for a flat circular 
washer made of steel (G = 0.8x10° MPa, v = 0.25, p= 7.8 
kg/m*) with a 2:1 ratio of outside radius to inside radius, 
its inner contour under fluctuating uniformly distributed 
pressure and the outer contour free. Calculations made 
for three frequencies of pressure fluctuations (1250 s°', 
2500 s', 5000 s'') and for three values of the damping 
coefficient a (0.002pw, 0.02pw, 0.2pw) indicate that an 
increase of the pressure fluctuation frequency results in a 
higher initial energy, while an increase of the damping 
coefficient, too, resuits in a higher initial energy but also 
in a longer decay time and a lower final energy. Figures 
2; references |. 


Action of Acoustic Wave on Set of Two Parallel 
Circular Cylinders in Ideal Fluid 

917F0164G Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 1, Jan 91 pp 117-123 


[Article by A.P. Zhuk, Institute of Mechanics, UkSSR 
Academy of Sciences, Kiev] 


UDC 534.2:532 


[Abstract] Action of an acoustic wave on two parallel free 
solid circular cylinders, one behind the other in the path 
of the wave, 1s analyzed by treating it as a problem of 
diffraction. The problem is formulated for a plane 
acoustic wave and two such cylinders with generally 
different radii r, , in a boundless medium, a compress- 
ible ideal fluid. The problem is formulated in a rectan- 
gular global system of coordinates Oxyz with the origin 
at a point half-way between the parallel axes of the two 
cylinders, the Ox axis in the direction of wave incidence 
crossing them at right angles, and the Oz axis running 
half-way between and parallel to them. The problem 
reduces to calculation of the velocity potential in the 
wave field, a mathematically linear problem of solving 
the linear wave equation, and to calculation of the 
radiation force. Solving the problem of diffraction for 
the incident primary wave with a given velocity potential 
requires calculating the velocity potential of the sec- 
ondary wave field, which must satisfy the conditions of 
radiation at infinity and boundary conditions compat- 
ible with those for the primary field at the surfaces of 
both cylinders. This part of the problem 1s formulated in 
two cylindrical local systems of coordinates, one for each 
cylinder with the origin on its axis. It is solved, after 
separation of variables, in series of Hankel functions. 
Owing to symmetry with respect to the xOz plane, both 
cylinders oscillate with their velocity vectors in the Ox 
direction. The problem now reduces to an infinite series 
of linear algebraic equations solvable only by the trun- 
cation method. The radiation force on the cylinders 1s 
calculated by time-averaging the hydrodynamic force of 
the fluid medium on them. This force and the velocity 
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potential must be calculated with the same accuracy, 
including quadratic terms in parameters of the wave 
field. The action of an acoustic wave with a given 
amplitude (0.9 cm?/s or given intensity (175.5 W/m* and 
frequency (8 kHz) on two cylinders in propanol (density 
0, = 0.7854 g/cm’) and the dependence of their attendant 
interaction on the distance between them are analyzed 
numerically according to this method. First are consid- 
ered two cylinders with different radii (r, = 1.3r,) in two 
configurations: smaller cylinder behind the larger and 
larger cylinder behind the smaller. The dependence of 
the radiation force at each cylinder on that distance 
between them 1s found to follow a slightly different trend 
in each case, evidently owing to different distributions of 
the secondary wave field. A larger radiation force is 
tound to be always acting on the larger cylinder and the 
radiation force acting on the second cylinder, larger or 
smaller, to be equal to the radiation force acting on 1t if 
there were no other cylinder before it. The first cylinder 
thus does not influence the second one, which 1s consis- 
tent with the fact that reflected acoustic energy propa- 
gates in the opposite direction. Cylinders with equal 
radii are found to approach each other, faster or slower 
depending on the magnitude of the parameter kr (k- 
wave number. Calculations for this special case were 
made assuming an L/r = 4 length-to-radius ratio. The 
difference between the radiation forces acting on the two 
cylinders and causing them to approach each other does, 
moreover, becomes larger with higher frequency of the 
incident acoustic wave. Figures 5; references 5. 


Resonance Dynamics of Vessel With Liquid 
Mounted or Vibration Dampers 


917F0164F Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 1, Jan 91 pp 109-117 


[Article by G.N. Puchka and V.V. Kholopova, Institute 
of Mechanics, UkSSR Academy of Sciences, Kiev] 


UDC 629.78 


[Abstract] The response of a cylindrical vessel with 
liquid mounted on stationary elastic vibration dampers 
to a periodic external excitation 1s analyzed for the 
important practical case of motion in a field of linear 
potential forces, 1.e., where nonlinear effects of inertia 
dominate over nonlinear effects of elasticity. The equa- 
tions which describe the resulting oscillations of the 
vessel and the free liquid surface are derived from the 
Hamilton-Ostrogradskiy principle of least action and 
formulated in two Cartesian systems of coordinates, one 
of them “absolutely” stationary (O*X*Y*Z*) and one 
rigidly affixed to the vessel (OX axis vertical and coin- 
ciding with c/linder axis). The height of the liquid 
column ts assumed to be smaller than the height of the 
cylinder (vessel not full) but equal to or larger than its 
radius. For a determination of the velocity potential in 
the liquid, the boundary-value problem for the appli- 
cable equations of motion can then be solved by Nari- 
manov's method (G.S. Narimanov, PRIKLADNAYA 
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MATEMATIKA I MEKHANIKA Vol 75, 21, No 4). The 
vessel is free to move in both orthogonal horizontal OY 
and OZ directions. As to oscillations of the free liquid 
surface, the first two antisymmetric modes (f,, and f,) are 
most significant in terms of energy content and only they 
are being considered. The system of equations of motion 
is resolved with respect to higher-order derivatives, 
third-order derivatives of motion parameters having 
been retained in it, for subsequent analysis of the vessel- 
liquid dynamics by the method of averaging. A change 
from the natural space coordinates Y,Z to quasi-normal 
ones Y = X,+ X>, P = ¥,X,+ Yox> and Z = X4+ X4, F = ¥4X, 
+ Y4x, facilitates determination of the four natural 
frequencies @, > in the two directions in the the XOY 
plane and @, , in the two directions in the XOZ plane. 
The system of equations of motion is resolved with 
respect to higher-order derivatives, the third-order deriv- 
atives of the motion parameters having been retained, 
for subsequent analysis of the vessel-liquid dynamics by 
the method of averaging. The conditions are thus estab- 
lished under which the two antisymmetric oscillation 
modes of the free liquid surface will appear and when 
oscillations of the vessel-liquid system in one plane only, 
the XOY plane of excitation, will change to oscillations 
in both XOY and XOZ planes. Special attention is paid 
to such compound oscillations in two orthogonal planes 
in the special case of an attendant internal resonance, 
with two pairs of equal natural frequencies w, = w, and 
®> = W,, in addition to resonance excitation at a fre- 
quency ® = @, = @, OF @, = @,. Figures 3; references 6. 


Stability Theory for Motion of Triple-Unit 
Vehicles on Tires 


917F0164E Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 1, Jan 91 pp 96-103 


[Article by L.G. Lobas, V.P. Sakhno, and T.I. Tarnopol- 
skaya, Institute of Mechanics, UkSSR Academy of Sci- 
ences, Kiev, and Kiev Institute of Automobile Roads] 


UDC 531.384:629.114 


[Abstract] A mathematical model describing plane- 
parallel motion of triple-unit motor vehicles on tires 1s 
constructed for a stability analysis of that motion. As a 
specific example is considered a tandem trailer truck 
consisting of a tractor and two trailers, the first one being 
a platform without wheels which clears the road surface 
Change from Cartesian space coordinates x.y of the 
inertial system to quasi-velocity coordinates v = (dx/ 
dt)cosé + (dy/dt)sind, u = -(x/dt)sind + (dy/dt)cosd (8- 
course angle of tractor) makes it possible, by virtue of the 
laws kinematics, to separate the variables in the equa- 
tions of dynamics and thus split the system of these 
equations into two lower-order subsystems of two equa- 
tions for forces and moments. In addition to the driving 
and reaction moments there are also included here 
damping and restoring moments which viscoelastic sup- 
porting-and-coupling devices produce, these moments 
depending on the trailer folding angles o, , and their 
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time derivatives. The resulting linear combination of 
four equations is reformulated in variations, for a sta- 
bility analysis of uniform rectilinear motion describable 
by an inertia matrix in accordance with the single-track 
model of a tandem trailer truck which in this case 
includes two supporting- and-coupling devices. An ana- 
lytical solution for the critical velocity is supplemented 
with numerical data pertaining to typical! such trailer 
truck. The critical velocity is shown to depend on the 
distance from the tractor-platform coupling point to the 
center of mass of the tractor, becoming infinitely high 
when that distance 1s the critical one and becoming zero 
when that distance is smaller. Figures 4; references 12. 


Vibrations of Reflector Antenna Shells Under 
Vibratory Excitation 

¥17F0164D Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 1, Jan 91 pp 64-71 


[Article by V.S. Gudramovich, N.G. Baranov, N.A. 
Konovalov, and I.M. Maytala, Institute of Engineering 
Mechanics, UkSSR Academy of Sciences, Dnepro- 
petrovsk] 


UDC 538:539.2:773.5.004 


[Abstract] An experimental study of reflector shells for 
large antennas was made concerning their vibrations 
under vibratory loads, such loads acting during seismic 
activity On ground antennas as well as on spaceborne 
antennas while the spacecraft 1s put in orbit by a rocket 
A paraboloidal 0.35 mm thick shell serving as the larger 
reflector of a two-reflector antenna, made of nickel 
(modulus of elasticity E = 7x10* MPa, tensile strength 
140 MPa, 0.2 % yield point 175 MPa), was tested in both 
vertical and honzontal positions on a VEDS-1 500 elec- 
trodynamic exciter. The length of the shell was 0.44 m 
and and the radius of its base was 0.78 m. The shell was 
tested in a reinforcing hoop and without it, the cross- 
section of the hoop being 0.0115 m° and also without 
large. For a determination of its amplitude-frequency 
characteristics and vibration resistance, it was tested 
over the 5-4000 Hz frequency range under 6-12 g accel- 
erations from 0.6 g to 12 g. The vibration amplitudes 
were measured with five ANS-14 piezoelectric acceler- 
ometers spaced 0.180 m apart on a flange made of 
Al-Mg(6M) aluminum alloy and mounted around the 
shell at the 0.231 m radius of the latter. The frequency 
marker, developed at the Institute of Engineering 
Mechanics, was a device consisting of two pairs of 
stationary light-emitting photodiodes and a rotating disk 
with two rows of holes coupled to an excitation fre- 
quency indicator. The electric output signals from those 
accelerators were recorded on an NO-67 magnetograph 
and also on an NO-45 loop oscillograph, with the AGK-2 
harmonic analyzer developed at the Institute of Eng- 
neering Mechanics as matching device. The vibration 
modes of the shell with and without hoop at two reso- 
nance frequencies, 9 Hz and 140 Hz, were photographed 
with a “Gladiolus-!"’ camera at the maximum rate of 


Mechanics of Gases, Liquids, Solids 15 


240 frames/s with a 90° shutter so as to allow 1/190 s 
exposure time, and recorded on KN-4 film. The results 
of this study indicate, most importantly, the effective- 
ness of stiffening such an antenna reflector shell. Figures 
5: tables |; references 4. 


Design of Viscoelastic Thin-Walled Shells 
Considering Their Damageability 


917F0164C Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 1, Jan 91 pp 55-61 


iArticle by Ya.A Eyubov, Azerbaian Institute of Struc- 
tural Engineering, Baku] 


UDC 539.3 


[Abstract] A method of designing viscoelastic thin- 
walled shells with consideration of their damageability 
and geometrical nonlinearity is developed by application 
of the variational principle to calculation of stresses and 
Strains, this principle being derived from the one for a 
three-dimensional body. Solution of the corresponding 
boundary-value problem is difficult, because the con- 
straint on the solution is the differential equation of 
damage process kinetics and the functional depends on 
the damageability parameter. Since only a limited 
number of coordinate functions will satisfy this con- 
straining differential equation of kinetics, an exact solu- 
tion of it is waved by replacing it with an approximate 
integral relation. The functional, whose Euler equation 
contains the kinetics of damage, thus becomes one of the 
mixed kind. The nonlinear problem of viscoelasticity for 
given boundary condition and then for the initial condi- 
tions is solved accordingly, by the Ritz method. As an 
example is considered a cylindrical viscoelastic shell, a 
pipe, hinge-supported at both end and under an internal 
pressure uniformly distributed over the lateral surface. A 
numerical analysis of the solution to the corresponding 
Cauchy problem indicates that increasing the wall thick- 
ness of such a shell will lengthen the critical time tll 
fracture occurs, but not proportionally. Figures 2, refer- 
ences 4. 


Use of Acoustoelasticity Theory of Rayleigh 
Surface Waves for Determination of Stresses in 
Solid Bodies 

917F0164B Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 1, Jan 91 pp 44-49 


[Article by A.A. Chernoochenko, F.G. Makhort, and O.1 
Gushcha, Institute of Mechanics and Institute of Electric 
Welding imeni Ye.O. Paton, UkSSR Academy of Sci- 
ences, Kiev] 


UDC §$39.3.:620.179.16 


[Abstract] Biaxial stresses which electric welding produces 
in the surface layer of an elastic and initially isotropic solid 
body are calculated according to the acoustoelasticity 
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theory of Rayleigh surface waves in the linear approxima- 
tion. The sum and the difference of principal stresses o, , + 
6,, and 6,,- 6,, are first expressed in terms of the 
velocities of those waves propagating through that layer 
along the principal axes OX, and OX,. For a better 
interpretation of experimental data, they are then 
expressed in terms of the recirculation frequencies of 
ultrasonic waves or pulses used for measurements by the 
recently developed acoustic method. In an experiment 
with two disks, one of D16 aluminum alloy 35 mm thick 
with a 180 mm in diameter and one of 10Mn2VB alloy 
steel 15 mm thick with a 120 mm diameter, biaxial stresses 
in them were produced by diametral compression along 
the vertical OX, axis with a constant load of 100 kN 
(aluminum alloy) and 400 kN (steel) respectively. Stresses 
were measured along the horizontal diameter on the OX, 
axis at various distances from the center. Ultrasonic sur- 
face waves in these disks were transmitted and received by 
a pair of piezoelectric transducers, changes 1n the recircu- 
lation frequencies during the compression process being 
recorded by instrument in the 10~-error accuracy class. 
Each transducer consisted of two 4 mm thick and 10 mm 
wide plates made of TsS-19 (Sr zirconate-titanate) piezo- 
ceramic, rigidly joined to form a ““wedge”’ with a resonance 
frequency of 3 MHz. In a separate test rectangular plates 
20 mm thick, 50 mm wide, and 70 mm long were com- 
pressed lengthwise while the recirculation frequencies of 
ultrasonic surface waves were measured along and across 
the plates. The results validate this method of stress 
measurement. It therefore was used for measurement of 
residual stresses in an 8 mm thick, 70 mm wide, and 310 
mm long rectangular plate hard-faced by electric welding. 
for calibration of the proportionality factors in the theo- 
retical formulas. Figures 3; tables |; references |! 


Stressed and Strained State of Half-Space Due to 
Translatory Motion and Rotation of Partly 
Embedded Rigid Ball Partly 


9¥17F0164A Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 1, Jan 91 pp 24-31 


[Article by I.K. Senchenkov, V.G. Savchenko, O.P 
Chervinko, and V1. Bobyr, Institute of Mechanics, 
UkSSR Academy of Sciences, Kiev] 


UDC $39.3.01 


[Abstract] A stress and strain analysis 1s performed for a 
composite material which consists of an elastic binder 
matrix with fine abrasive filler grains moving in it, and 
the integral stiffness characteristics of such a system are 
evaluated in connection with the use of such materials 
for cutting tools. The model for this analysis is a per- 
fectly hard spherical grain only partly embedded in an 
only slightly deformable isotropic and linearly elastic 
matrix occupying a half-space. The grain rotates about 
its two axes respectively parallel and perpendicular to 
the matrix surface, while it also moves in a straight line 
but so constrained that its embedment depth remains 
constant. Interaction with neighboring grains 1s assumed 
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to be negligible. The depth of grain embedment in the 
matrix 1s characterized by the ratio A = A’/R, (R,- radius 
of grain), 0 < A < 2 (d= 0 corresponding to a grain on 
the matrix surface and A = 2 corresponding to a grain 
completely under the matrix surface). The motion of a 
spherical grain is treated as that of a solid body under a 
force N =(N,.,N,,N,) and a moment M = M, ,M,.M.,). In 
view of the axial symmetry of the problem, a system of 
cylindrical coordinates is added for displacements. The 
equations of balance for the components of force and 
momentum are then solved for boundary conditions of 
zero displacements on the matrix surface and zero 
Stresses on the grain surface. T'ie nonzero boundary 
conditions for strains are stipulated on points of the 
circular corner along which the grain surface intersects 
the matrix surface. The problem is solved by the method 
of finite elements, its solution being sought in the form of 
a trigonometric series with respect to the angular circum- 
ferential coordinate. A numerical analysis of the solution 
reveals how not only the stresses and the strains in the 
matrix but also of the elements of its stiffness matrix 
depend on the grain embedment depth, including their 
asymptotic behavior as the that depth tends to infinity. 
Figures 5; references 5. 


Stressed State of Noncircular Thick-Walled 
Cylinders Depending on Indicators of Load 
Nonuniformity Along Generatrix and Directrix 


917F02174 Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 3, Mar 91 pp 12-18 


[Article by Ya.M. Grigorenko, G.G. Viaykov, and S.N 
Shevchenko, Institute of Mechanics, UkSSR Academy of 
Sciences, Kiev] 


UDC $39.3 


[Abstract] The stressed state of noncircular thick-walled 
hollow cylinders 1s analyzed for its dependence on the 
degree of load nonuniformity along the generatrix and the 
directrix. The analysis is premised on the three- 
dimensional theory of elasticity. A noncircular cylinder 1s 
placed in an orthogonal curvilinear system of coordinates 
s, t, y (coordinates s = const and t = const on lines of 
principal curvatures of the median surface, representing 
families of generatrices and directrices respectively, coor- 
dinate y in direction normal to the median surface). The 
problem 1s formulated as a system of six partial differential 
equations 6V/t = F(V.8V/5s,6V/5y.6°V/bs° 8° V 

ady’ 8°V/Ss5y,.P.T) (P- load vector, T-temperature field) 
which has been resolved for the first derivatives of its 
components with respect to the t-coordinate. Satisfying the 
conditions of periodicity or symmetry of the stressed state 
will reduce solution of the problem for closed cylinders to 
solution of the problem for open ones, these conditions to 
be satisfied by a compatible combination of sufficiently 
smooth functions defining the boundary conditions at 
sections of a cylinder bounded by its directrix t = const 
Boundary conditions at lateral surfaces are stipulated in 
terms of stresses. The problem 1s solved in the region 0 < 
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s<L.t<t<ty,-H/2 < y <= H/2(H, =1,H,=R,,h, 
= |, R,- radius of curvature of directrix which defines the 
median surface). The problem is solved by combining 
separation of variables with subsequent direct and discrete 
orthogonalization. First the dimei.sionality of that system 
of equations is lowered by expansion of given and sought 
functions into Fourier series in the s-coordinate. Separa- 
tion of variables then reduces the three-dimensional 
problem to a series of two-dimensional ones for each 
Fourier harmonic. This series of two-dimensional prob- 
lems in s,y coordinates is then, by the straight lines 
method, reduced to a one-dimensional boundary-value 
problem for an N*-order system of ordinary differential 
equations (N* = 6n- 4, n- number of points on the 
finite-differences approximation over the cylinder wall 
thickness) dV*/dt = A(t)V*+ B(t) where A, V*(t, = B, and 
A,V*(t,. = B, (A, 5- rectangular N*/2 X N* matrices, B(t)- 
given N*-dimensional vector, B,, - given N*/ 
2-dimensional vectors). This equation can now be solved 
very accurately by the numerical method of discrete 
orthogonalization. This 1s demonstrated on hollow ellip- 
tical thick-walled cylinders. The perimeter of the median 
surface 1s equated to the circumference of a circle with a 
radius R and the cross-section 1s described by two para- 
metric equations x = b.cos@, z = a.sin® (a,b- semiaxes of the 
ellipse). The results indicate an interaction of loads non- 
uniformly distributed along the two coordinates s and t 
respectively. Figures 4; tables 2; references 5. 


Solution of Problems Pertaining to Parametric 
Vibrations of Shells of Revolution by Method of 
Finite Elements 

Y17F0217B Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 3, Mar 91 pp 32-37 


[Article by A.T. Vasilenko, S.S. Kokoshin, and P.N. 
Cherinko, Institute of Mechanics, UkSSR Academy of 
Sciences, Kies 


L Dc $39.3 


[Abstract] Parametric vibrations of thin-walled shells of 
revolution with a uniform thickness are analyzed by the 
method of finite elements, the invariance property of 
this method making it possible to solve the problem of 
such vibrations for shells with essentially any geomet- 
rical and mechanical characteristics under any con- 
straints and under any mechanical or thermal loads. As a 
specific example ts considered a shell under a harmonic 
perturbation force. The functional of total potential 
energy for this problem consists of four terms, each a 
surface integral, two of them representing respectively 
the generalized strain energy in the perturbed state at the 
instant of stability loss and the strain energy during 
variation of displacements under constant stresses in the 
subcritical state. A third term representing kinetic energy 
is added and a fourth term representing work done by 
the load being subtracted. The relation between general- 
ized strains € and generalized displacements U 1s 
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expressed in the form e¢ = LU (L- matrix of differentia- 
tion operators). Following separation of variables for the 
displacement vector, the kinematic relations between 
Strains € and displacements of nodal points A, are for 
each finite element are expressed in the form € = BA, (B- 
matrix which includes the L matrix). Angular displace- 
ments w are related to linear displacements U through a 
matrix L’ of differentiation operators so that w = L’U 
and to displacements of nodal points A, through a matrix 
G (which includes the L’ matrix) so that m = GA,. The 
expression for the potential energy of each finite element 
is then expressed in terms of A, displacements rather 
than strains. The coefficients with which these displace- 
ments appear are the physical stiffness matrix, the geo- 
metrical stiffness matrix, the mass matrix, and the vector 
of generalized time-dependent actions (forces, moments) 
at nodal points. All three matrices, each a double inte- 
gral, are evaluated numerically by the Gauss method. 
The condition of steadiness of these functionals with 
respect to A,” (velocity) and A, “ (acceleration) leads to 
two systems of algebraic resolvent equations det] ...| = 0 
and thus to an eigenvalue problem for two systems of 
equations. In the absence of external loads and for 
correspondingly zero initial stresses, moreover, the 
problem reduces to natural vibrations of a shell. Other- 
wise the problem is characterized by involving two 
parameters, namely initial stresses as load parameter in 
addition to frequency. It can be solved by inverse itera- 
tions with the aid of Rayleigh’s relations. The method 
was applied to two specific problems, one a test problem 
with a known analytical solution and one which had 
already had been solved numerically by the method of 
discrete orthogonalization. The first problem was to 
calculate, in accordance with the classical theory of 
shells, the spectrum of lower-order natural vibrations of 
a closed cylindrical one freely supported at both ends. 
The second problem was to analyze parametric vibra- 
tions of an isotropic cylindrical shell fixed at one end in 
a sliding constraint at the other. Tables 2. references | 2. 


Axisymmetric Problem of Impact of Solid Body 
against Thin Membrane Floating on Compressible 
Fluid Half-Space 

917F0217C Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 3, Mar 91 pp 38-45 


[Article by V.D. Kubenko, Institute of Mechanics, 
UkSSR Academy of Sciences, and V.V. Gavrilenko, 
Institute of Automobile Roads, Kiev] 


UDC 533.6.013.42 


[Abstract] Impact of an axisymmetric blunt solid body of 
otherwise arbitrary shape against a thin circular mem- 
brane with a finite radius floating on the boundary of a 
compressible fluid which occupies a half-space 1s consid- 
ered, the body having moved in the direction of its axis 
of symmetry and approached the membrane normally to 
its surface before striking it at the center. The fluid 1s 
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assumed to be a weightless barotropic ideally compress- 
ible one. The problem is formulated in the linear approx- 
imation, with corresponding constraints on the velocity 
v,(t) of the body and on the deflection w* of the median 
surface at its center after impact at time t = 0: v,(t) <v 
(speed of sound in the fluid) and w* < R (radius of 
membrane). The shape of the body surface at any instant 
after impact is described by the equation z = ®t,8) (6- 
polar angle from axis of body and membrane symmetry 
to membrane-fluid contact cucle, ®- slightly curving 
piecewise-smooth even function). The wave function for 
the fluid and the classical equations of dynamics for the 
membrane are modified so that their solutions can be 
obtained in the form of series in the same system of 
Legendre polynomials, while the boundary conditions 
are linearized and then reduced to those at the quiescent 
fluid surface. This is made possible by identifying the 
linear coordinates along the body, membrane, and fluid 
surfaces for the early stage of body-...embrane-fluid 
interaction. In order to solve the problem then, it 1s 
necessary to first solve two auxiliary ones. One of them is 
establishing the relation between the strain rate of the 
fluid surface V(t,8) and the hydrodynamic pressure p(t,8) 
on it, which is done by Laplace transformation of both 
and subsequent expansion of their Laplace transforms 
into a Fourier series in Legendre polynomials. The other 
problem ts establishing the relation between both deflec- 
tion w(t,8) and rate of deflection dw(t,6)/dt and the 
pressure p(t,6), which 1s done analogously and by 
applying the convolution theorem for two original func- 
tions. On the basis of this theoretical model, the problem 
of impact was solved numerically for a mgid sphere and 
a thin steel membrane on water. Figures 4; references 8. 


Viscoelastic State of Cylindrical Triple-Layer 
Shell Under Pulsed Heat or Force Load 

91 7F0217D Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 3, Mar 91 pp 46-52 


[Article by E.I. Starovoytov and S.A. Vorobyev, Belorus- 
sian Institute of Railroad Transportation Engineers] 


UDC $39.3 


[Abstract] A not so flat cylindrical shell with a carner- 
filler-carner structure and an asymmetncally nonuni- 
form thickness under a pulsed heat or force load 1s 
considered, the material of the filler being different than 
the material of the two carners but both materials being 
isotropic. Both matenals are charactenzed by elasticity 
of volume changes and a linear viscoelastic lag in shear 
The behavior of the two outer layers conforms to the 
Kirchhoff-Love hypotheses and that of the inner one 
conforms to the exact theory of elasticity in the linear 
approximation of displacements from the transverse 
coordinate. The boundary conditions at the contact 
surfaces are continuity of displacements. Strains are 
assumed to be small. The problem is formulated in a 
cylindncal system of coordinates and solved for tangen- 
tial displacements and normal ones (deflections) of 
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points on the median surface of each layer. It is solved by 
application of the Hamiulton-Ostrogradskiy principle. 
which leads to a system of six integrodifferential equa- 
tions describing the motion of such a shell. This system 
of equations is solved by the Bubnov-Galerkin method, 
first for the quasi-static problem as a test case. In this 
case it reduces to a system of Volterra equations of the 
second kind, after the inertia terms have been discarded 
on account of their smallness. The solution of this 
problem establishes the confidence level for analysis of 
the shell dynamics, that problem being solved by te 
numerical Denisenko-Yanovich method (N.V. 
Denisenko and L.A. Yanovich, DIFFERENTSIAL- 
NYYE URAVNENIYA Vol 19 No 5, 1983) for 
boundary conditions of free supports at both ends. The 
method is based on reversibility of the differential oper- 
ator associated with three constant diagonal! matrices yp, 
x, p. As a first example is considered a normally incident 
uniform thermal flux heating the shell continuously from 
time t = 0 on. The transient temperature field nonuni- 
formly distributed over the stack thickness is then calcu- 
lated according to Starovoytov’s formula, disregarding 
the curvature of layers. As a second example is consid- 
ered heating of the shell for a certain time interval only 
and subsequent stabilization of the temperature field. 
Creep is included in numerical calculations for carner 
layers of DI6T duralumin and a filler layer of teflon 
(polytetrafluoroethylene). As a third example 1s consid- 
ered an impulse A&t) of pressure (& Diurac’s delta- 
function) loading the shell axisymmetrically. Figures 4; 
tables |; references 9. 


Elastic Displacements Vector of Thin Spherical 
Shell 

91 7FO217E Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 3, Mar 91 pp 63-69 


[Article by V.A. Merkulov, Volgograd Institute of Con- 
struction Engineering} 


UDC $39.3:534.] 


[Abstract] Considering that the displacement vector of a 
closed thin spherical shell needs to be known for a 
precise analysis of its natural vibrations, its elastic dis- 
placements vector U and elastic rotations vector I of the 
median surface are calculated on the basis of a linear 
displacement profile across the shell wall thickness in 
accordance with the Kirchhoff-Love hypotheses. This 
profile is described in the form of the vector relation U, 
= LU + 2zF xe, (U,- displacements vector, z - thickness 
coordinate, x - vector product. e, - normal unit vector) 
To this relation are added relations describing the linear 
dependence of strain components on the distance from 
the median surface, these relations corresponding to 
Love's equations of state. Differentiation of the general- 
ized displacement vector with respect to each of the two 
geographical coordinates generates a system of two dif- 
ferential equations relating the two derivatives to the 
radius of the median surface of the shell in terms of 
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strains and angular displacements for the vibration 
mode vector V,. Equality of its mixed derivatives being 
the condition for integrability of that system of equa- 
tions, they together with the relations describing linear 
strain profiles across the shell wall thickness yield rela- 
tions between angular displacements and strains. These 
relations are, in accordance with Love’s equations of 
State, transformed into a system of three equations of 
balance of forces (normal forces T,+ T,= T, shearing 
force S) and moments (bending moments M,+ M,= M, 
twisting moment H) which in turn yields a system of 
three resolvent equations for the three u,v,w displace- 
ments and then a system of two equations for natural 
frequencies: one for the frequencies of vibrations of an 
element of the median surface about the normal to it and 
one for two series of frequencies of conjugate radial- 
tangential vibrations. Both sernes of frequencies w,, are 
then calculated by Novichkov’s method (Yu.N. Nov- 
ichkov, Handbook: “Vibrations in Engineering”, Izd. 
Masshinostroyeniye, 1978). One senes is [q<]infn 
with n = 0,1,2,.... the n = 0 mode being the tension- 
compression mode of equal purely radial displacements 
of points of the median surface. The highest frequencies 
of tangential tensile and shearing vibrations have been 
calculated according to the Novichkov’'s formula with A, 
+ (A.- 4B_'? in the numerator for this series. The other 
series 1S @, with n = |, 2,... and no n = 0 mode. The 
lowest frequencies of radial tensile and shearing vibra- 
tions have been calculated according to Novichkov’s 
formula with A,- (A,- 4N_)'* in the numerator for this 
series. References | | 


Numerical Analysis of Statics of Smooth and 
Discretely Reinforced Shells Under Combined 
Flexural and Axial Load 

YI17F0217F Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 3, Mar 91 pp 69-75 


[Article by V.P. Maksimenko and N.V. Kovalchuk, 
Institute of Mechanics, UkSSR Academy of Sciences, 
Kiev] 


UDC §39 3 


[Abstract] The static stressed and strained state of closed 
thin-walled circular cylincrical shells discretely rein- 
forced by an orthogonal grid of hoops and stringers, 
stiffer rings at both ends, under a combined axial and 
flexural load 1s analyzed following a solution of both 
Viasov’s equation of balance and compatibility equation 
for strains (V.Z. Vlasov, “General Theory of Shells and 
Its Engineering Application”, Goztekhizdat, 1949). The 
load at the upper end of the shell is an axial force q 
uniformly distributed over the base circle and a concen- 
trated radial transverse force Q applied to the top ring 
The load at the lower end of the shell is four axial forces 
P distributed each over one of the four s = nR/6 long arcs 
(R- radius of shell) symmetrically spaced around the base 
circle. These are the reactions of four mounts on which 
the shell rests, spaced symmetrically 90° apart under the 


Mechanics of Gases, Liquids, Solids 19 


lower base either so that one pair lies on the diameter 
parallel to the line of action of force Q and the other pair 
lies on the diameter in quadrature or so that both pairs 
lie on diameters at 45° to the line of action of force Q. 
The problem is for each case split into two separate ones, 
for the transverse load and for the axial load respec- 
tively, this being permissible by virtue of the principle of 
independent action of forces. The first problem has one 
plane of symmetry and one plane of antisymmetry. The 
second problem has only several planes of symmetry. 
Both problems are also solved for smooth shells, as a 
basis for evaluating the effectiveness of reinforcement. 
The equations of balance and compatibility, each with 
four finite series on the left-hand side, are reduced to a 
system of linear algebraic equations for numerical solu- 
tion by the Gauss method. As an example is considered 
a shell made of Al-Mg(6M) aluminum alloy (radius R = 
20 cm, length L = 5xR/8, wal! thickness t = R/400) with 
32 stringers, 4 hoops, and two end rings. The results are 
compared with experimental data, tests having been 
performed on 42 such shells with a different Q/P load 
ratio for each. Figures 5; tables 3; references 7. 


Supercritical Deformation of Real Cylindrical 
Shells Under External Pressure and Estimation of 
Their Stability 

917F0217G Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 3, Mar 91 pp 76-83 


[Article by A.Yu. Yevkin and V.L. Krasovskiy, Dnepro- 
petrovsk Institute of Construction Engineering} 


UDC $39.3 


[Abstract] Formulas for the sensitivity of shells to per- 
turbations based an analysis of supercritical deformation 
according to Coyter’s [transliteration] theory were on 
checked out experimentally on real cylindrical shells 
under external pressure for a verification of stability 
estimates. Special shells for this experiment were pro- 
duced from strips of normalized Cr!18N19 stainless steel 
by spot welding them into cylinders with an R = 71.5 mm 
inside radius, an active length L = 2.0R, and a wall 
thickness t = 0.00375R. Each shell was hermetically 
closed at both ends so as to prevent leaks into the cavity 
and ensure constancy of the constraints, which approxi- 
mated those of free support. External loading was 
effected by sucking air out of the cavity. In the first 
experiment preliminary depressions simulating local 
imperfections were produced either so as to leave no 
residual stresses or so as to conform to the stable 
supercritical configuration with one local depression. In 
the second experiment preliminary protrusions simu- 
lating local imperfections were produced by buckling 
under various under various lower pressures. For a 
second evaluation of the effect of subcritical regular 
imperfections “resonating” with the supercritical shell 
configurations on the magnitude of the critical pressure, 
the results of this experiment are interpreted by consid- 
ering the equilibrium position of a shell after large 
deformation. Theoretically this problem 1s 
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solved on the basis of nonlinear equations in Volmir’s 
mixed theory of shells (A.S. Volmir, “Stability of 
Deformat ie Systems”, Izd. Nauka, 1967). It is solved by 
the asymptotic method using Pade’s approximations. A 
comparison of stability estimates according to this and 
Coyter’s method indicates that the latter are erroneous in 
the case of large deformation. This method therefore is 
preferable when initial deflection must be taken into 
account, which can be done by expressing the depen- 
dence of the load parameter on the total deflection in the 
form a complete quadratic polynomial. Figures 5; refer- 
ences 1 2. 


Vibrations of Triple-Layer Plate of Finite 
Dimensions in Contact With Liquid 


917F0217H Kiev PRIKLADNAYA MEKHANIKA 
in Russian No 3, Mar 91 pp 88-93 


[Article by S.N. Beshenkov and T.D. Volkova, Zapor- 
ozhye University] 


UDC 539.3:534.1 


[Abstract] A triple-layer plate of finite dimensions fas- 
tened to an an infinitely stiff frame and floating on the 
surface of a liquid which occupies a half-space is consid- 
ered, its forced harmonic vibrations being analyzed in 
two ways. First they are analyzed on the basis of the 
simplified theory for a sheath-filler-sheath plate which 
disregards transverse compressibility of the filler in 
accordance with the broken-line hypothesis. Its deforma- 
tion is described by the Goloskokov-Filippov system of 
three equations (Ye.G. Goloskokov and A.P. Filippov, 
“Nonsteady Vibrations of Deformable Systems”, Izd. 
Naukova Dumka, 1977) in deflection w of the plate, 
tangential displacement u of the median plane of the 
filler, and angular displacement w of its plane sections 
under a net load equal to the difference q - p between 
excitation force and reaction of the liquid, the time 
factor e‘’' being omitted. This system of equations is 
solved for boundary conditions of free support at the 
ends with w, u, w as well as gq and p expressed in the form 
of infinite Fourier series. The longitudinal profile of the 
acoustic pressure, reaction of the liquid, is approximated 
as a piecewise-linear one and its unknown magnitude of 
that pressure at the nodes is calculated using the Huy- 
gens integral of plate deflections. Insertion of w, u, and 
into that integral then leads to a system of algebraic 
equations for reaction of the liquid at the nodes. The 
problem is also solved on the basis of the refined theory, 
taking into account motion of the filler and describing it 
in accordance with the dynamic theory of elasticity: c,’ 
grad div r - c,* curl curl r + w’r = 0 (c, - velocities of 
longitudinal and transverse waves in filler, r- displace- 
ment vector, w- vibration frequency). Normal and 
shearing stresses in the filler along its boundaries with 
the two sheaths are calculated in accordance with 
Hooke’s law. The continuity condition for displacements 
in the plate then yields the deflections of the lower sheath 
in contact with the liquid. The acoustic pressure of the 
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liquid is calculated using the Huygens integral of these 
deflections, which again leads to a system of algebraic 
equations for reaction of the liquid at the nodes. Both 
methods were applied to two plates with a 2h = 2 mm 
thick filler of a material with a modulus of elasticity E = 
5x10’ N/m? in one plate between two aluminum sheaths 
of the same thickness t = 3 mm (symmetric structure) 
and in one plate between two aluminum sheaths of 
unequal thicknesses t, = 5 mm, t, = | mm (asymmetric 
structure), both plates being excited by a linear force F = 
1000t(x- L/2) N/m (L- length of plate). Calculations were 
made for a -umerical evaluation of the dependence of 
the plate deflection at its center and of the reaction of the 
liquid (acoustic pressure) at that point on the excitation 
frequency. Both deflection and pressure were treated as 
complex quantities, their real and imaginary parts being 
calculated separately. Deflection of the symmetric plate 
structure and acoustic pressure on it were found to have 
analogous frequency characteristics, with maximum 
imaginary part and zero-crossover of the real part at the 
same resonance frequencies (150 s' and 1410 s). The 
results of calculations for the asymmetric plate structure 
reveal the effect of that asymmetry on the two frequency 
characteristics. The results of calculations by the second 
method, moreover, reveal the effect of compressibility of 
the filler. Figures 3; references 8. 


Interaction of Acoustic Wave and Two Parallel 
Cylinders in Viscous Fluid 

9]7F02171 Kiev PRIKLADNAYA MEKHANIKA 

in Russian No 3, Mar 91 pp 108-115 


[Article by A.P. Zhuk, Institute of Mechanics, UkSSR 
Academy of Sciences, Kiev] 


UDC 5$32.516:539.376 


[Abstract] Interaction of an acoustic wave and two 
parallel free circular cylinders with generally different 
radiir, >, one behind the other in the path of the wave, in 
a viscous compressible fluid is treated as a problem of 
diffraction. The problem is formulated in a rectangular 
global system of coordinates for the velocity potential of 
the wave field and in two cylindrical local systems of 
coordinates, one for each cylinder with the origin on its 
axis for the two successive diffractions of the incident 
primary acoustic wave. The rectangular system of coor- 
dinates, its origin at a point O half-way between the 
parallel axes of the two cylinders, is onented with the Ox 
axis in the direction of wave incidence crossing them at 
right angles and the Oz axis running half-way between 
and parallel to them. The mean force of the fluid on each 
cylinder is defined as a surface integral, the product of 
the stress tensor in the fluid by the normal unit vector 
being integrated over the surface of each cylinder, and is 
calculated by time-averaging this integral. The mean 
force acting along the Ox axis was calculated numerically 
for two cylinders with equal radii R = 0.1 cm in propanol 
(density p, = 0.7845 g/cm’), L = | cm apart, and an 
acoustic wave of 175.5 W/m? intensity. The results 
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indicate that the mean force on each cylinder depends on 
the wave frequency and on the density p of its material. 
When its density 1s lower than that of the fluid and the 
cylinders float, then the first cylinder influences the 
behavior of the second more than the second cylinder 
influences the behavior of the first. In both cases of p, > 
p and p, < p the two mean forces acting on each cylinder 
(force of the incident wave and force of the wave 
reflected by the other cylinder) act in opposite direc- 
tions, the two cylinders attracting each other in the first 
case aud repelling each other in the second case. As the 
viscosity of the fluid decreases with rising temperature, 
moreover, an acoustic wave produces a smaller mean 
force on the cylinders. Figures 5; references 6. 


Design of Framed Shells of Intricate Shapes for 
Internal Pressure Load 
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{Abstract} An economical algonthm is proposed for 
design calculation of the stressed and strained state of 
linear shells in reinforcing frames, geometrically linear 
structures, taking into account deflections of the joints 
with the reinforcing members. It 1s based on the method 
of finite elements using the Hess matrix H" of second 
derivatives of the potential energy with respect to gener- 
alized coordinates, triangular finite elements are selected 
for a thin shell. Their generalized elastic displacements 
vector 1s defined in terms of the strain tensor of their 
median surfaces and the generalized q-coordinates 
vector q 1s defined in terms both linear and angular 
displacements. For finite elements of the reinforcement 
members both vectors are defined analogously, but refer- 
ring to corners of a finite element and disregarding the 
nonlinear terms. The generalized q-coordinates q 1s 
determined from the system of equations H''g = R (R- 
vector of external forces). The algorithm has been pro- 
grammed for a BESM-6 high-speed computer. It was 
tested on shells of noncanonical shapes conforming to 
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the Kirchhoff-Love hypotheses under an internal pres- 
sure load and the results compared with experimental 
data, the agreement being quite satisfactory. Figures 5; 
references 7. 


Experimental Stability Study of Nonconservative 
Systems 
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[Abstract] The test stand built for experimental stability 
study of nonconservative systems where the pressure of 
the air jet discharged by a nozzle at the free end of a 
uniform cantilever beam serves as the tracking compres- 
sion force P on this beam was used for determining the 
critical length L of the beam according to the P., = 
14.8EJ_,;,, /L’ criterion rather than the P., = 20.2EJ,,,./ 
L? based on the dynamic criterion of stability. The 
tracking force in this test stand is produced by a pro- 
peller driven by a small variable-speed 0-600 rpm 1.47 
kW electric d.c. aircraft motor weighing 2.6 kg. The 
tested beam, a rectangular bar of quenched 25CrMnSiN, 
steel, 1s mgidly clamped at one end to the test stand table 
in a horizontal cantilever position. The motor with an 
aerodynamically streamlined propeller on the shaft 
extension is suspended on a thread so as to face the free 
end of the beam rather than being attached to it. This 
arrangement eliminates deflection of the beam. The 
twisting moment, which the propeller produces in addi- 
tion to the axial tracking compression force, is absorbed 
by a carriage also mounted at this end of the beam and 
consisting of a two slotted guide bars for a steel ball each 
inside a yoke. The tracking force can be varied by 
varying the voltage at the motor terminals. The beam 
can thus be subjected to a force from 0 to 20 N 
continuously and up to 30 N for a short time. One 
drawback of the test stand are small vibrations, due to 
some rotor and propeller unbalance, these vibrations 
becoming larger and their effect on the test results 
becoming more appreciable with increasing motor 
speed. Tests were performed during the 1986-89 period 
by many researchers including M.Kh. Mullagulov, on 
four beams 70,80,90, 100 cm long with a rectangular 0.97 
cm x 3.94 cm cross- section. The beams were also tested 
for stability under a tracking tension force, attempts to 
induce instability in the form of either whip or flutter not 
having been successful. Figures 2; tables |; references 5. 
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Design of Closed-Loop for Control of Robot With 
Torque Motors 
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in Russian No 12, Dec 90 pp 15-20 


[Article by A.V. Timofeyev, V.A. Baykalov, and I.R. Fom- 
ina, Leningrad Institute of Aircraft Instrument Design] 


UDC 62.50 


[Abstract] Automatic control of the manipulator in 
robots with torque motors is synthesized on the basis of 


the Lagrangevector differential equation describing the 
manipulator dynamics, proportional-differential regula- 
tors being considered for this system. Stability of regu- 
lator operation has been selected as the criterion for 
design and selection of the drive for each degree of 
freedom. Accuracy and response speed of the control 
system are then estimated for design review. The calcu- 
lations have been programmed on a Standard System 
1045 computer in various languages (FORTRAN, 
PASCAL, PL/1) for educational and research purposes. 
A typical menu is shown for >omputer-aided design of 
such an automatic control for SKARA robots. Tables |: 
references 5. 
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